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Abstract 

We derive functional equations for the eigenvalues of the XXZ model subject to anti-diagonal 
twisted boundary conditions by means of fusion of transfer matrices and by Sklyanin's method of 
separation of variables. Our findings coincide with those obtained using Baxter's method and are 
compared to the recent solution of Galleas. As an application we study the finite size scaling of 
the ground-state energy of the model in the critical regime. 
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I. INTRODUCTION 

The Quantum Inverse Scattering Method (QISM) [l^ provides a powerful framework for 
the construction of exactly solvable lattice models with the Yang-Baxter equation defining 
the underlying algebraic structure. This structure also serves as the basis for the solution of 
the spectral problem, i.e. explicit construction of eigenvalues and eigenstates, for the transfer 
matrices of these lattice models, by means of the algebraic Bethe ansatz. The application of 
this method, however, is limited to systems with a simple pseudo vacuum state which can be 
used as a reference state to generate the complete spectrum by application of the elements 
of the Yang-Baxter algebra. There exist methods to compute the spectrum of these systems 
which do not suffer these restrictions, e.g. Baxter's method of commuting transfer matrices 
{2] or Sklyanin's method for separation of variables [3;]. In these approaches, the eigenvalues 
are encoded into the solutions to certain functional relations which have to be solved in a 
second step: using the analytical properties of the transfer matrix its eigenvalues can be 
parametrized in terms of the roots of Bethe equations. 

Recently, such alternative routes to the solution of the spectral problem of exactly solvable 
models have been applied to anisotropic spin chains with open ends and general non-diagonal 
boundary fields. Here, the boundary fields break particle number conservation and therefore 
the completely polarized state of spins cannot be used as a pseudo vacuum for the algebraic 
Bethe ansatz. Various methods have been used to determine the spectrum of the correspond- 
ing transfer matrix 0, [sl, sl, [?, I^, I^, 10|. Some of these methods require restrictions on the 



boundary fields and/or the bulk anisotropy of the spin chain to work, other approaches rely 
on conjectures regarding the existence of certain limits. Depending on the specific approach 
one obtains very different functional equations for the spectral problem and in many cases 
an efficient procedure for their solution is still missing. 

In this paper we consider some of these approaches in a simpler setting to better under- 
stand how they are related and to assess their applicability for the solution of the spectral 
problem: the XXZ model with general toroidal boundary conditions is defined by the hamil- 
tonian 

L 

^ = E H^^i + ^l^^i + coshr/a|a|^.i] , a^^^ = R-'a^K (1.1) 
where a", a = x,y,z denote the Pauli matrices for spins \ at site j. The unitary matrix K G 
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End(C^) determines the boundary conditions. For anti-diagonal K the model is integrable 
but has no pseudo vacuum state (see below) . It has first been solved by means of Baxter's 
method and solutions to the resulting functional equations can be given in terms of the 
roots of Bethe ansatz equations. 

The integrability of (11. ip is established within the QISM [l] from the Yang-Baxter algebra 

/?i2(A - ^)Ti(A)r2(^) = T2(/x)Ti(A)i?i2(A - /i) . (1.2) 

Here the monodromy matrix Tj(A) is a matrix on the auxiliary linear space Vj with entries 
being the generators of the quadratic algebra. The structure constants are arranged in the 
i?-matrix, Rjk{X) G End(V^- V^) which itself solves the Yang-Baxter equation 

i?12(A - fx) i?i3(A - U) i?23(/i -l^)= i?23(/i ' 1^) Riz{\ - v) RuiX - /i) • (1-3) 

In this paper we will use the well-known trigonometric solution for two-dimensional spaces 
Vj corresponding to the six-vertex model. In this case a representation of (11. 2p is given by 
local ^-matrices with Tq = Raj = Cj reading 

^^^^^ /sinh(A)^o + cosh(A)^| Sr \ (1.4) 

^ \^ 5+ sinh(A)50 - cosh(A)5|y ' 

The elements of Cj are operators on a two-dimensional quantum state space of a spin |, 
in terms of the Pauli matrices aj they are S'O = cosh(r//2)/, = sinh(r//2)(T| and Sf = 
sinh(?7/2) cosh(r7/2)crj'^. Using the co-multiplication of the Yang-Baxter algebra new repre- 
sentations can be constructed from (II. 4p : the monodromy matrix T(A) = Cl{\) ■ ■ ■ Ci{\) 
also satisfies (II. 2p . Note that Cj{X) (and T(A)) satisfies 

CjiX + ivr) = -</:,(A)ao^ = -a]Cj{X)a] (1.5) 

where the subscript denotes the auxiliary space. 

Apart from these operator valued representations there exist 2x2 c-number solutions K 
to (II. 2p . Being independent of the spectral parameter A they satisfy [-R(A), K(^K~\ =0. For 
the six-vertex i?-matrix, this relation has two classes of solutions, namely diagonal or anti- 
diagonal twist matrices K. Without loss of generality {12!. [l3| we restrict our considerations 
to 

K=\ \ , K=\ \ (1.6) 





with a twist angle in the diagonal case. 

As a consequence of fll.2l) the transfer matrix t{X) = tYQ[KT{X)] generates a family of 
commuting operators, [t{X),t{fi)] = 0. Therefore the spin chain hamiltonian (11. ip which is 
obtained from 

d\nt{X) 



H = 2 sinh r] 



A-2 
^-2 



— Lcoshrj (1.7) 



dX 

is integrable. 

In the following section we will first briefly review the existing solutions of the XXZ model 
with anti-diagonal twisted boundary conditions followed by two different approaches to the 
spectral problem based on (A) the fusion hierarchy of transfer matrices [14!] at anisotropies 
rj = m/{p + 1) with integer p > 1 and (B) Sklyanin's method of separation of variables j^, 
respectively. In Section UTI] we study the spectrum of small chains to identify the ground-state 
and extract the critical properties from the finite size scaling behaviour of the ground-state 
energy in the massless regime rj = i'^ with real 7. 

II. SOLUTION OF THE SPECTRAL PROBLEM FOR ANTI-DIAGONAL TWIST 

For a diagonal twist matrix the eigenvalues and eigenstates of the transfer matrix t{X) 
can be obtained by means of the algebraic Bethe ansatz starting from the ferromagnetic 
so-called pseudo vacuum with polarized spins, (cr|) = 1 (see e.g. [1|). For the anti-diagonal 
twist matrix 

however, the total magnetization is not a conserved quantum number. As a consequence 
there is no simple reference state such as the ferromagnetic one and the algebraic Bethe 
ansatz cannot be applied. Instead, a functional equation (called TQ-equation) for the eigen- 



values A(A) of the trans 



transfer matrices 
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12| 



'er matrix has been obtained using Baxter's method of commuting 



A(A)g(A) = sinh^ (a + |) g(A - r^) - sinh^ (a - |) g(A + r/) . (2.2) 
This difference equation is solved by 

g(A) = nsinh-(A-A,). (2.3) 



As a consequence of the analyticity of the transfer matrix eigenvalues it follows that the 
rapidities Xj are different solutions to the Bethe equations 

sinli^(A + i,,) „ smhl(A-At + .;) f^^A^ai' (2 4) 

smh-(A-l,)--ll.inhl(A-A.-„)' '"^^(H^i' P'"' 

(Note that these equations with an extra phase (—1) and any number M < L of rapidities Xj 
determine the spectrum of a staggered six- vertex model 15|]. This case, however, cannot be 
obtained from (11. ip with the twist matrices (II. 6p .) From (ll.7p we obtain the corresponding 
eigenvalue of the spin chain hamiltonian (11.11) 

/ r . s sr^ sinh ri sinh ^ 

^»^^>'=^-"" + ^S: eo.hA,-coshr P-^) 

J = l 

Recently, Galleas has proposed a different approach to solve the spectral problem of 
the spin chain with anti-diagonal twist. From the Yang-Baxter algebra he derives a closed 
set of equations for the L — 1 zeroes A[.^^ of the eigenvalues A(A) and a second set of 2(L — 1) 

(2) 

rapidities A) defined to be the zeroes of the matrix element of a certain element of the 
Yang-Baxter algebra between the ferromagnetic state and the eigenstate of the model with 
twist 



sinh(A« + i,) sinh(Ai^) - i,) ^^ ' ' _ e^ ^^^^^^ sinh(Af - X^^) 

l\ sinh(Af-ir/) 

sinh(A^^^ - A^ 
t\ sinh(A« - Iv) 




L-l 



2i7rr 



(2.6) 



Here r = 0, . . . , 2L — 1 is a quantum number fixing the translational properties of the 
corresponding eigenstate. Note that as a consequence of the functional TQ-equation (12. 2p 
A[.^^ are the 'hole solutions' A = A^,^'' ^ {Xj}^^^ to the Bethe equations (12.40 . 

The algebraic equations (12. 6p involving two sets of rapidities are vaguely reminiscent of 
so-called nested Bethe equations for systems with higher rank symmetry. In the present 
case, however, the two sets of rapidities are redundant in the sense that - for an eigenstate 
with given r - either set can be eliminated in favour of the other one. Also, unlike the string 
hypothesis for the Bethe rapidities Xj nothing is known about the structure of the solutions 
of Galleas' equations in the thermodynamic limit L ^ oo, which would be a prerequisite for 
both numerical or analytical studies of large systems. 
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A. Fusion hierarchy and truncation identity at roots of unity 



This method was first developed for the RSOS model by Bazhanov et al. [IJ] and was 



adapted to spin chains by Nepomechie e.g. la, ll7|]. Unfortunately this method only works 



if the values of the crossing parameter are chosen to be roots of unity r] = i7r/(p + 1) 



Nevertheless as in the periodic case 17|] the solution obtained is valid for arbitrary t] as it 
coincides with (12. 2p . 

Because of the Yang-Baxter algebra higher spin transfer matrices of the XXZ spin chain 
obey a so-called fusion hierarchy, i.e. it is possible to construct a higher spin transfer matrix 
out of lower spin transfer matrices directly. A transfer matrix using a spin-j auxiliary space 
is constructed using fused ^-matrices. 

The fused spin-(j, i) ^-matrix for j = i, 1, |, . . . is given by |19| 



^^■" = ;^E^^ ■ (2-8) 



/:<i...2,>2,+i(A) = P+..2,£i,2,+i(A)£2,2,+i(A + v)--- A„2,+i(A + (2j - l)r7)P+..2, , (2.7) 

where is the projector defined by the sum over all permutation operators for n spin-| 
spaces 

(7 

In the same way the fusion of the twist matrix (12.11) is carried out yielding the spin-j 
representation 

ir<i...2,> = P^U,KiK2 ■ ■ ■ K,^Pt.2, ■ (2.9) 

The fused monodromy matrix for the chain then reads for L lattice sites of the original 
hamiltonian (II. ip 

T<i...2,>(A) = /:<i...2,>,l(A) ■ ■ ■ /:<i...2,>,i(A) (2.10) 
and tracing out the auxiliary space gives the associated transfer matrix 

tO)(A) = tri...2,i^<i...2,>T<i...2,>(A) . (2.11) 



For these transfer matrices the following fusion hierarchy with j = |, 1, |, . . . holds 17|, ll9|. 
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2l| 



tO+V2)(^) = tO)(A)t(i/2)(A + 2j7]) - idgT){X + (2j - l)7])t^^-'^^\\) , (2.12) 
where t^^/^\X) = t{X) and t^^\X) = I. 



On 
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the other hand fused transfer matrices can be constructed using quantum-group theory 



22| . The £- matrices of higher auxihary spins from quantum-group constructions have a 



simple direct relation to an ^-matrix with lower auxiliary spin at roots of unity rj = i7r/(j9+l) 
with p being an integer number. It is also possible to relate the quantum-group /^-matrices 
to those constructed via fusion. Resulting in the identity at roots of unity 



Bi...p+iAi...p^iC^i...p^i^^p+2WAi}.p-i-iBilp^i 



/i(A) 



\ 



(2.13) 



where the entries of matrix A are unnormalized Clebsch- Gordon coefficients in the decom- 
position of the tensor product of 2j spin-1/2 representations into a direct sum of SU(2) 
irreducible representations and the matrix i? is a diagonal matrix needed for symmetrizing 

Q. 

The function yu(A) is related to the quantum determinant 



KT)(A-r/) = -(-/i(A))^ 
/i(A) = — sinh(A + ^rj) sinh(A — ^rj) 



(2.14) 
(2.15) 



and z/(A) = — /i(A)~^(i/2)P sinh ((p + 1)(A — ^i])) is related to the crossing parameter via p. 

The fused twist matrices themselves obey a truncation identity similar to fl2.13p . Under 
the fusion procedure an anti-diagonal matrix with only I's as entries remains anti-diagonal 
with I's as entries after applying the transformation of the appropriate Clebsch-Gordon 
matrix and omitting null rows and columns, hence 



A,...,,K. 



2i-f^<l-2j> 



All 



2i 



^l--2j-2-^<l-2j-2>^1.^.2j-2 



A 



(2.16) 



Identities (12.131) and (I2.16P together give a truncation identity for the product of twist 
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and monodromy matrix 



Bi...p+iAi...p+iKT^i...p^iy^p+2i^)^i}.p+iBilp^i 



( 



i?i...p_iAi...p_iKT<i...p_i> p(A + r])AY,,,p_^B^,,,p_^ 



L 



(-KA))"F^ 



(2.17) 



with F = Ylj=i and accordingly for the transfer matrix by taking the trace of f l2.17p 

t(p+i)/2(A) = -{dgT){X - 7]){-lf&-^^/\X + T]) . (2.18) 

The fusion hierarchy f l2.12p together with the truncation identity (12.181) leads to a func- 
tional relation for the transfer matrix at roots of unity for a given p, e.g. for p = 2 or j = 1 
respectively this relation is 

t(A)t(A + r])t{X + 27]) - {d^T){X)t{X + 2r/) + (rf,T)(A + ri)t{X) + 

+ {-l)\d,T){X-r])t{X + ri) = ^. (2.19) 

Like in the RSOS model [l^ or the periodic XXZ chain [itI the goal is to recast the 
general form of the functional relation (I2.19P as a determinant of a certain matrix. This 
determinant being zero ensures the existence of a null eigenvector which leads to equations 
similar to TQ-equations. 

In the case of an anti-diagonal K-matrix the functional relation found above cannot be 
recast directly, though multiplying it with itself shifted by ivr = {p+l)ri results in a recastable 
expression. For general p this is a determinant of a (2p + 2) x [2p + 2) matrix reading with 
the eigenvalue A of the transfer matrix t 



( 



det 



Ao 





ho 
Ai hi 
-h As 



-irh,\ 





—h2p+i Asp 











V(-i)^H+i 

In the above expression we used the shorthands 

Afe =A(A + kri) 
/ijt = sinh^(A + kr] 



(2.20) 



-/iq Asp+i / 



2' 



(2.21) 
(2.22) 
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The definition of directly reveals = (— l)'^/ip+i+fc. This and the periodicity of the 
eigenvalue = — (— l)^Ap+i+fc following from (11.51) are needed to verify the equivalence of 
the determinant and the product of functional relations. 

Let (go, • • • , '?2p+i) be the null eigenvector of the matrix, this yields the equations 

Aogo + hoqi - {-l)^hiq2p+i = 

-hk+iQk-i + Afc^fe + huqk+i = for A; = 1, . . . , 2p (2.23) 
(~l)'^^2p+i9o " hoq2p + A2p+ig2p+i = . 

Using the ansatz = q{X + ki]) with 

^ 1 

g(A) = nsmh-(A-A,) (2.24) 

i=i 

the equations (12.231) imply only a single TQ-equation 

A(A)g(A) = sinh^(A + Ir]) g(A - r/) - sinh^(A - |r/) g(A + r]) (2.25) 

agreeing with (12.21) and leading to the same Bethe ansatz equations (12. 4p . Notice the 27ri 
periodicity of the g-function arising from the 2(p + 1) rows of the matrix in (I2.20p and the 
product in (I2.24p running up to L due to the structure of the upper right and lower left 
entries. 



B. Separation of variables 

In this section, we carry out the procedure of separation of variables, generalizing 
Sklyanin's result for the XXX chain 

We modify our definition of the monodromy matrix by introducing inhomogeneities 5j 

T(A) = CLiX -SL)---CiiX- Si) . (2.26) 

Later we will discuss the limit 6j 0. We employ the usual notation for the elements of 
the twisted monodromy matrix 



A{X) B{X) 

ax) D{x) 



KT{X) . (2.27) 
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The Yang-Baxter algebra contains the commutation relation 

[B{X),B{f,)]^0. (2.28) 

It is therefore reasonable to assume that there exists a complete set of A-independent eigen- 
vectors \i) of B{X). Note that this does not follow from the commutativity since we have not 
shown yet that B{X) is diagonalizable. By induction in L we show that in the standard basis 
D{X) is lower triangular with all diagonal entries equal to zero, and B{X) is lower triangular 
with diagonal entries 

sinh(A -xi)--- sinh(A - 4) where = 5^- ± f . (2.29) 

This shows that -B(A) is indeed diagonalizable, provided the inhomogeneities Sj are mutually 
distinct, and that its eigenvalues are given by 

B{X) \e) = sinh(A - xi) ■ ■ ■ sinh(A - x[) {£) . (2.30) 

In other words, we have defined operator-valued zeroes Xj of -B(A) 

B{X) = sinh(A - xi) ■ ■ ■ sinh(A - xl) , (2.31) 

where Xj = diag(a;], . . . ,x'j^) in the eigenbasis of B{X). Since the joint spectrum of the 
operators Xj is not degenerate, any eigenvector \£) is completely determined by its eigenvalues 
, . . . , x^. We interpret the set of eigenvalues (xf, . . . , x^) of a given eigenvector \i) as a 
point in C^. Then the Hilbert space of the spin chain is isomorphic to the space of complex- 
valued functions on the set X C of these points (for this reason the separation of variables 
method is also called 'functional Bethe ansatz'). In this picture the Xj are the operators of 
multiplication by the coordinate functions Xj in C"^: 

xjf = x,f , (%/) (xO = x'^fix') (2.32) 

for any function f on X. In the following, we shall not distinguish between the operators 
Xj and the functions Xj. 

We want to formulate the spectral problem for the twisted transfer matrix in the diagonal 
basis of the operators Xj. To this end, we first define the 'conjugated momenta' X^ to the 
'coordinates' Xj, which we obtain from ^(A) and D{X) by substituting A = Xj 'from the left' 

{i\X-\m) = {i\A{x'j)\m) and {e\X+\m) = {i\D{x'j)\m) . (2.33) 
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The following commutation relations hold, which can be shown in the same way as for the 
XXX case 



[xj,Xk] = , 
XpXk = (xk ± ri^jk)Xf 

[X+,X,1=0 for j^k, 
XfXj = A{x, ± f ) , 



(2.34) 



(2.35) 



(2.36) 



(2.37) 



(2.38) 



where A (A) is the quantum determinant of the twisted monodromy matrix 




These commutation relations largely fix the action of the conjugated momenta Xf" on the 



eigenvectors of the operators xj. The remaining freedom is due to the fact that the xj- 
eigenvectors are determined only up to phase factors [sl, Thm. 3.4]. (In the functional 
language changing these phase factors corresponds to multiplying all functions with an 
arbitrary function which has no zeroes.) We will now reconstruct the conjugated momenta 
from the commutation relations. This will allow us to formulate the spectral problem for the 
twisted transfer matrix in the Xj-eigenbasis without knowing the change of bases explicitly. 
Let uj be the function lu = 1 on X, and define the functions by 




(x) = {Xfu;){x) . 



(2.40) 



Equation f l2.35p implies 




(2.41) 




(xx, . . . , Xj, ■ ■ ■ y XXj^ 



{xi,...,Xj±r],...,XL) . 



(2.42) 



As the functions / are defined only on X, 




whenever 




x^X 



(2.43) 
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must hold. The commutation relations fl2.36l) - (l2.38l) translate into the following conditions 
on the functions A^: 

Af{x)A^{Epx) = A^{x)Af{Eix) , (2.44) 
A+(x)A-(E;x)=A-(x)A+(E-x) for j^k, (2.45) 
A±(x)Aj(Efx)=A(x,±f). (2.46) 

We make the following ansatz: Let the functions A^ be defined as 

Af(x) = A±(x,) , (2.47) 

where 

A±(A) = e± sinh(A - 5i ^ f) ■ ■ ■ sinh(A -5lT^) ■ (2.48) 

In this definition, the constants ^± are an arbitrary factorization of the determinant of the 
twist matrix 

e+e- = det(K) = -1 , (2.49) 

and the functions Aj. factorize the quantum determinant of the monodromy matrix in the 
following sense: 

A+(5,-f)A_(5, + f) = A(5,). (2.50) 

This ansatz satisfies the conditions given in Eqs. fl2.43p - fl2.46p . 

We now return to the spectral problem for the twisted transfer matrix 

t{X)q{x) = A{\)q{x) . (2.51) 

The eigenvector q is in our representation a complex-valued function on X. We substitute 
A = Xj from the left and obtain 

Aixj)q{x) = A+{xj)q{Ejx) + A_{xj)q{E-x) . (2.52) 



The coefficients A± now depend on only one coordinate Xj. This is due to the particular 
solution Af to the conditions ^T^ - ^IM^ chosen in Eqs. ^TM^-^IM^, for a generic solution 
it would not have been the case. The last equation (12. 520 suggests the separation of variables 
ansatz 

q{x) = qiixL) ■ ■ -qiixi) . (2.53) 
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It remains to solve a set of one-dimensional problems 

A{xj)qj{xj) = A^{xj)qj{xj + r])+A^{xj)qj{xj-r]) , j = l,...,L, (2.54) 

which we recognize as the TQ-equation fl2.2p . evaluated on the discrete lattice X. 

Recalling that Xj takes the values 6j ^i]/2 and using Aj-{6j ±i]/2) = 0, we see that the 
finite-difference equation fl2.54p takes the form of a homogeneous system of linear equations 



A(5, + f) -A_(5, + f)\ /g,(5, + f) 
-A45,-f) A(5,-f) ) ^(5,-f) 

For a nontrivial solution its determinant has to vanish 



. (2.55) 



A(5, + f )A(5, - f ) = A{5,) , J = 1, . . . , L (2.56) 

(we have used equation fl2.50p ). The eigenvalue A (A) is of the form 

A(A) = A„i+ie(-^+^)^ + A_L+3e(-''+=^)^ + ■ ■ ■ + Az^.ie^^-^)^ . (2.57) 

The L coefficients Aj have to be determined from the L equations (12.561) . Each of these 
equations defines a quadratic form in the L-dimensional complex space of the coefficients of 



A (A) (see also lOj). These quadratic forms intersect at 2^ points, which correspond to the 
2^ eigenvalues A(A). (Note that the eigenvalues of the transfer matrix are non-degenerate 
for the anti-diagonal twist. This explains why the method of separation of variables does 
not suffer the so-called 'completeness problem' of the algebraic Bethe ansatz. In the case 
of diagonal twist matrices the method cannot be applied since the operator -B(A) has only 
L — 1 zeroes.) 

Finally, we want to remove the inhomogeneities Sj, which were not present in the original 
problem related to the spin chain fll.ip . Simply putting them all to zero is not possible since 
then the equations f l2.56p are no longer independent. Instead, we consider the equation 

A{6 + f )A(5 - f ) = - sinh^((5 - r/) sinh^((5 + r]) (2.58) 

and its derivatives w.r.t. 6 up to order L — 1, evaluated at 6 = 0. For small systems we have 
verified that this procedure does give the correct eigenvalues A (A). 
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III. FINITE-SIZE SCALING OF THE GROUND-STATE ENERGY 



A. Solution of the Bethe equations 

As an application of the solution of the spectral problem we now study the description of 
the ground-state in terms of the solutions to the Bethe equations. We start by solving these 
equations (12 ■4p for the complex Bethe rapidities {Aj}jL^ numerically for small systems. 

In the antiferromagnetic massive regime, described by real r], the Bethe rapidities of the 
ground-states, and only those, are purely imaginary for small system sizes and small values of 
1]. As the system size grows or for larger values of rj two of the rapidities in the ground-state 
form a single '2-string' of rapidities symmetric to the imaginary axis, A ± 7^/2. Computing 
the energy of this state one obtains a system size independent contribution. This is the 
energy of a domain wall, corresponding to the interfacial tension in the six-vertex model 
computed in Ref. [ll 



Here we concentrate on the massless case, i.e. rj = i'-f with real 7. In Figures [T] and 
[2] some results for up to 4 lattice sites and rj = itt/A are shown. As a consequence of the 
periodicity of the Bethe equations the set {Xj + ivrj^^i is also a solution which parametrizes 
a second state corresponding to a different eigenvalue (12. 2p of the transfer matrix but has 
the same energy fl2.5l) . As the lattice size increases, we find that all Bethe rapidities of the 
two ground-states have imaginary parts or tt. Below we will use this observation to study 
the finite size scaling behaviour of the ground-state energy. For excited states we have not 
identified such a pattern which would be necessary for a systematic analysis of the excitation 
spectrum starting from the Bethe equations (12.41) . 

Also shown in Figures [T] and [2] are the corresponding solutions to Galleas' equations (12.61) . 
i.e. the two sets of rapidities {X^j^^}j~i, {Xf'^}'^j[!^^^^ and the quantum number r. Galleas' 
equations are invariant under separate shifts of the two sets by ivr, therefore each rapidity 
appears twice and A and A -|- ivr have to be identified. In this parametrization, the pairs of 
degenerate solutions are given by the same rapidities A^-^'* and X^^ , but their values of r differ 
by L. Again, the distribution of the roots of (12.61) in the complex plane does not appear 
to follow a simple scheme as the system size increases, not even for the the ground-state. 
Therefore, this approach is of limited use only to study the spectrum of long chains. 

To proceed, we parametrize Bethe roots corresponding to the ground-states as follows: 
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FIG. 1: Distribution of the Bethe rapidities Xj (•) and the rapidities A^^^j (□) and A'-^^j (v) in the 
massless regime with rj = in /A for two lattice sites. 

let /ij, j = 1, . . . , m be the real Bethe rapidities and uj, j = 1, . . . , n the real parts of the 
Bethe rapidities with imaginary part it. For even system sizes L we have m = n = L/2, for 
odd system sizes we have m = {L^ l)/2, n = (L± l)/2 for two ground-state configurations. 
We take the logarithm of the Bethe equations (12. 4p and obtain 



fc=l k=l 
m n 



m 



(3.1) 



e=i 



where 



$(A) = -iln 
e(A) = -iln 
e(A) = -iln 



sinh(A - ^) 
sinh(A + f ) 
sinh[l(A — i7)] 
sinh[i(A + i7)] 
sinh[i(A + ivr — 17)] 
Lsinh[i(A + i7r + i7)] 



TT + 2 arctan[tanh(A) cot(^)] , 



TT + 2 arctan[tanh(|) cot(|)] 



2/J ' 



2 arctan[tanh(|) cot(^ + |)] 



We have chosen the branch of the logarithm in such a way that the functions $, O and G are 
monotonically increasing. Here arctan denotes the principal branch of the inverse tangent 
taking values between —it/ 2 and it/ 2. Our numerical solutions for small chains show that 
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(a) L = 3, r] = m/A 



(b) L = 4, 7? = i7r/4 



FIG. 2: Distribution of the Bethe rapidities Xj (•) and the rapidities X^^^j (□) and A'-^^j (v) for 
the ground-state in the massless regime with r] = iir/A for (a) three lattice sites and (b) four lattice 
sites. 

the Bethe integers Ij and J^. in (13. ip follow a simple pattern: for even system sizes, one 
ground-state is obtained with 



while for odd L 



Ij=j-1 forj = l, 
Jk = k for k = 1, 

Ij = j - 1 for j = 1, . . 

for k = 1, . . 



L 

■'2 
L 

■'2 

L + 1 
2 

L - 1 



(3.2) 



(3.3) 



Jk = k 

The second ground-state is obtained by interchanging the sets {Ij} and {Jk}- 

In the thermodynamic limit we can solve the logarithmic Bethe equations (13.11) analyti- 
cally. Rewriting them as 

Yi^^,) = /, , Ziu,) = J, (3.4) 

with the 'counting functions' Y and Z which satisfy ¥{—00) = Z(— 00) = and Y{+oo) = 
m, Z{+oo) = n. Assuming that the distributions of the rapidities /ij and Vj can be described 
by countinuous densities p and a in the thermodynamic limit they are given by the derivatives 
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of the counting functions 



p(A) = Ir'(A) , a(A) = Iz'(A) . (3.5) 



(3.6) 



The logarithmic Bethe equations become a pair of coupled integral equations 
27rp(A) = ^'(A) - J dX' e'(A - A')p(A') - J dX' 0'(A - X')a{X') , 
2na{X) = <I>'(A) - J dX' e'(A - A')p(A') - J dX' e'(A - A')a(A') 

which can be solved by Fourier transform, yielding 

P(A) = ^(A) = ^ . . . (3.7) 

/7 cosh(A7r/7) 

Here we have assumed without loss of generality < 7 < vr. For the ground-state energy 
per lattice site in the thermodynamic limit we find with (12.51) 

sin7 f dX 2smj 
£oo = cos7 / — — — — . 3.8 

7 J cosh(A7r/7j cosh(zA) — COS7 

This agrees with the result for the untwisted chain |23|] , though the distribution of the Bethe 
rapidities is different. 

B. Conformal field theory 

The low energy effective field theory for the XXZ model in the massless regime, —1 < 
cos 7 < 1, is well known to be that of a free boson with compactification radius ^J2txR = 



^/(vr — 7)/7r. Conformal field theory predicts the finite size scaling of the energies of the 

25l | for a system with periodic boundary 



ground-state and the low lying excitations 
conditions as 

E,n{L)=Le^-^c+^{h + h)+ o{L-') (3.9) 

where Coo is the energy density in the ground-state, which is in our case given by Eq. (13.81) . 
and vf = ^ sin 7 is the 'Fermi' velocity of elementary excitations in the system. The 
universal number c is the central charge of the conformal field theory, for the free boson it 
is c = 1. 

From the energies appearing in the spectrum (13.91) for a given lattice realization of the 
field theory one can identify the operator content of the latter: h and h are the conformal 
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weights of primary operators of the CFT. By choosing particular boundary conditions one 
obtains different sectors of the theory (i.e. certain representations of the global symmetry 
group 0(2) of the model (11. ip ). For diagonal boundary conditions (11.61) with = and 7r/2 
the symmetry of the bulk Hamiltonian (II. ip is preserved and the spectrum is given in terms 
of the highest weights of two commuting U{1) Kac-Moody algebras. The scaling dimensions 
of the primary operators On.m are given in terms of the eigenvalues n of the U{1) charge 
operator 1 Ylj ci'^d momentum 2Txm/ L by 

2 

Xn,m = K,m + K,m = nRV + ^^ (3.10) 

(m takes integer (half odd integer) values for = (0 = 7r/2)). A diagonal twist (ll.6p 
with angle ^ 0,7r/2 breaks the global symmetry to 5*0(2). The dimensions of primary 
operators are Xn,m+^/n- As a consequence the finite size scaling of the lowest energy state 
energy is changed to Eqs{L) = Lcqo — t!'Vf Ccs/6L with an effective central charge [261] 

60^ 71 

Ceff = 1 - 12 Xo,^/^ = 1 z r, |0| < -• (3-11) 

TT^TT — 7j 2 

In the case of anti-diagonal twisted boundary conditions the 0(2) bulk symmetry is 
broken to Z2 (8> Z2 with the factors being generated by rotation around the 2-axis by tt and a 
global spin flip, respectively. The low-energy spectrum is that of a t/(l)-twisted Kac-Moody 



algebra without conserved charge. In this case the conformal weights are 



with integer fcj. 

We have solved the Bethe equations (13. ip for the ground-state of the spin chain numer- 
ically for system sizes up to L = 500. From (I3.12p we expect the lowest energy state to 
be that with conformal weights (/i, /i)o,o = {jq^ Jq) both even and odd number of lattice 
sites. With (13.90 this leads to the CFT prediction for the finite size scaling of the energy 

EosiL) = + ^ (4 + 12 (1^ + 1^) + . = 0- (313) 

The resulting effective central charge c^s = —1/2 is independent of the anisotropy which 
agrees nicely with our numerical data presented in Table [B In the isotropic limit of the XXZ 
model, 7 — i> 0, the spectrum depends only on the eigenvalues of the twist matrix (12.10 and 
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therefore the anti-diagonal twist is unitary equivalent to a diagonal one with twist-angle tt 

The corrections R{L) to the scaling (13.131) are a consequence of the fact that the lattice 
hamiltonian (11.11) differs from the conformally invariant hamiltonian of the continuum theory 
by terms involving irrelevant operators {27]. Perturbation of the conformal theory with an 
irrelevant operator with scaling dimension x > 2 leads to R{L) oc Therefore, by 

analyzing these corrections in the numerical data one can identify the leading irrelevant 
perturbation of the lattice hamiltonian. 

For the periodic XXZ chain and the model with diagonal twist the corrections to scaling 



vanish with an exponent x — 2 = max(47/(7r — 7), 2) [26|, [281]. The hamiltonian of the XXZ 
model with anti-periodic twisted boundary conditions, however, is related to the thermal 
operator of the Ashkin- Teller model and the leading corrections to scaling R{L) are generated 



by the operator Co,2 |29|, i.e. 



2-y 



R{L) oc L^-^°'2 = L"^ . (3.14) 

As shown in Figure El this provides an excellent fit for our numerical data. The dependence 
(I3.14P on the anisotropy parameter explains the slow convergence towards c^s in Table [T] as 
7-^0. 

IV. SUMMARY AND CONCLUSION 

For the XXZ spin chain with anti-diagonal twist which does not allow for a solution of the 
spectral problem by means of the algebraic Bethe ansatz due to the lack of a reference state 
we have derived the functional equations (12. 2p . originally obtained using Baxter's method 
of commuting transfer matrices, employing different methods: restricting the anisotropy to 



roots of unity rj = m/{p + 1) the TQ-equation of [ll| was obtained by truncation of the 
fusion hierarchy. In a second approach, namely through Sklyanin's separation of variables, 
we have used a representation of the Yang-Baxter algebra on a space of symmetric functions 
defined on a discrete lattice. In this formulation the spectral problem could again be recast 
in the form of the same TQ-equation (I2.54p - in this case, however, it has to be solved on 
the lattice of singular points of this functional equation only. This is unlike the situation 
for the spin chain with open boundaries and non-diagonal boundary fields where different 
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TABLE I: Finite size scaling of the ground-state energy in systems with even (odd) number 



of sites for several values of the anisotropy 7: Shown 


are the numerical results 


for —Ceff = 


{G/ttvf)L{Eo{L) 


— Lcoo) for lattices 


of size L together with the extrapolation to L = 


= 00. 


L 


7 = vr/S 


7 = 7r/4 


7 = 37r/8 


^ = tt/2 


30 


0.59536849 


0.53046883 


0.50367226 


500091 41 


60 


0.57644644 


0.51909229 


0.50156060 


50009285 


120 


0.56176871 


0.51200578 


0.50067152 


^0000^71 


240 


0.55014811 


0.50755893 


0.50029553 


U.(JUUUU1'4:0 


480 


0.54084022 


0.50476118 


0.50014993 


50000035 


extr. 


0.51(2) 


0.500(1) 


0.5000(1) 


500000 n") 


29 


0.40781510 


0.46933409 


0.49658570 


^000Q789 


59 


0.42572231 


0.48079424 


0.49850626 


50009363 


119 


0.43973453 


0.48794414 


0.49934855 


^00004Q4 


239 


0.45088912 


0.49242107 


0.49972106 


0.50000144 


479 


0.45987007 


0.49523121 


0.49990010 


0.50000036 


extr. 


0.48(2) 


0.499(1) 


0.4999(1) 


0.500000(1) 



functional equations have been found within different approaches. 

In the lattice formulation of the TQ-equation arising from the separation of variables the 
computation of the eigenvalues amounts to finding roots of L coupled polynomial equations 
f l2.56p . therefore this approach appears to be best suited to determine the spectrum of small 
systems. A similar limitation holds for the recent approach of Galleas j^. For an efficient 
solution of the TQ-equation in the thermodynamic limit the parametrization (12.31) or fl2.24p 
of the solution to the functional equations has to be used which leads to L algebraic Bethe 
equations (12. 4p . At least for the ground-state energy of the spin chain these equations can be 
solved, e.g. to obtain the interfacial tension of the six- vertex model (see [11]) or to identify 
the operator content of the low energy effective theory in the massless regime. 

Note, that together with the eigenvalues A(A) one obtains the functions q{X) from the 
TQ-equation in either approach. It is quite clear from the separation of variables that these 
functions contain the complete information on the eigenstates of the transfer matrix. How- 
ever, unlike the situation with the algebraic Bethe ansatz, where one has an expression for 
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10-^ 
10-^ 

^ -3 

— 10 



10-^ 



0.01 0.1 

1/L 

FIG. 3: Corrections to the scaling of the ground-state energy: Closed (open) symbols are numerical 
finite size data for even (odd) system sizes, lines are fits to a power law R{L) oc L~^'^/^'"~'^\ 

the eigenstates in terms of the generators of the Yang-Baxter algebra, the explicit transfor- 
mation from the g-functions to state vectors in the Hilbert space of the spin chain is not 
known. 
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